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ABSTRACT 


The  method  of  activation  energy  asymptotics  is  used  to  describe  the 
entire  history  of  a  carbon  particle  suddenly  immersed  in  a  hot  oxidizing 
ambient*  Under  appropriate  conditions  which  are  established  by  the  analysis 
the  history  of  the  particle  is  shown  to  consist  principally  of  a  heat**up 
period  during  which  no  significant  chemical  reaction  takes  place  followed  by  a 
period  involving  diffusion-limited  mass  loss  and  terminating  in  complete 
particle  consumption.  Between  these  two  extended  periods  brief  ignition  and 
post-ignition  periods  are  described#  The  final  demise  of  the  particle  occurs 
in  a  complex  extinction  period  during  which  all  the  quantities  describing  the 
particle  behavior  undergo  large  variations. 
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SIGNIFICANCE  AND  EXPLANATION 


Libby  and  Blake  treat  the  combustion  of  single  carbon  particles  in 
quiescent  cunbients  with  various  compositions.  Although  their  analysis  is 
applicable  to  a  variety  of  conditions,  those  encountered  in  entrained  flow 
combustion  are  emphasized.  Accordingly,  a  cold  particle  is  considered  to  be 
injected  into  a  hot  ambient,  generally  involving  oxygen  but  perhaps  other 
active  species  such  as  carbon  dioxide  and  water  as  well.  The  history  of  such 
a  particle  involves  an  initial,  pre-ignition  period  during  which  the  particle 
is  heated  but  undergoes  no  significant  chemical  reaction.  When  the  ambient 
temperature  is  sufficiently  high,  the  particle  reaches  a  critical,  ignition 
temperature  and  heterogeneous  reactions  involving  the  attack  of  carbon  become 
effective.  This  ignition  period  is  followed  by  post-ignition  behavior 
^.nvolving  rapid  increases  in  particle  temperature  which,  as  a  consequence  of 
chemical  reaction,  generally  exceeds  that  in  the  ambient.  Under  the 
circumstances  of  practical  interest  and  assumed  to  prevail  in  this  study  both 
the  ignition  and  post-ignition  periods  are  brief  but  result  in  the  rate  of 
mass  loss  by  the  particle  becoming  diffusion-limited.  There  ensues  an 
extended  period  involving  a  constant  rate  of  mass  loss  and  terminating  in  a 
complex,  although  brief  extinction  period  during  which  consumption  is 
complete.  The  brevity  of  the  transition  periods  associated  with  ignition  and 
extinction  implies  that  the  principal  features  of  particle  behavior  we 
consider  are  given  by  a  simplified  analysis  based  on  inert  particle  heat-up 
and  diffusion-limited  combusion.  The  present  study  employs  the  method  of 
activation  energy  asymptotics  (AEA)  to  rationalize  such  an  analysis,  to 
establish  the  conditions  for  its  applicability  and  to  resolve  the  structure  of 
the  transition  periods. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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INTRODUCTION 

Libby  and  Blake  (1-3]  treat  the  combustion  of  single  carbon  particles  in  quiescent  ambients  with  various 
compositions.  AJthougb  their  analysis  is  applicable  to  a  variety  of  conditions,  those  encountered  in 
entrained  flow  combustion  are  emphasized.  Accordingly,  a  cold  particle  is  considered  to  be  injected 
into  a  hot  ambient,  generally  involving  oxygen  but  perhaps  other  active  species  such  as  carbon  dioxide 
and  water  as  well.  The  history  of  such  a  particle  involves  an  initial,  pre-ignition  period  during  which 
the  particle  is  heated  but  undergoes  no  significant  chemical  reaction.  When  the  ambient  temperature  is 
sufficiently  high,  the  particle  reaches  a  critical,  ignition  temperature  and  heterogeneous  reactions  involv¬ 
ing  the  attack  of  carbon  become  effective.  This  ignition  period  is  followed  by  post-ignition  behavior 
involving  rapid  increases  in  particle  temperature  which,  as  a  consequence  of  chemical  reaction,  gen¬ 
erally  exceeds  that  in  the  ambient.  Under  the  circumstances  of  practical  interest  and  assumed  to  prevail 
in  this  study  both  the  ignition  and  post-ignition  periods  are  brief  but  result  in  the  rate  of  mass  loss  by 
the  particle  becoming  diffusion-limited.  There  ensues  an  extended  period  involving  a  constant  rate  of 
mass  loss  and  terminating  in  a  complex,  although  brief  extinction  period  during  which  consumption  is 
complete.  The  brevity  of  the  transition  periods  associated  with  ignition  and  extinction  implies  that  the 
principal  features  of  particle  behavior  we  consider  are  given  by  a  simplified  analysis  based  on  inert  parti¬ 
cle  heat-up  and  diffusion-limited  combustion.  The  present  study  employs  the  method  of  activation 
energy  asymptotics  (AEA)  to  rationalize  such  an  analysis,  to  establish  the  conditions  for  its  applicability 
and  to  resolve  the  structure  of  the  transition  periods. 
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The  method  of  A£A  is  a  recent  development  in  the  theory  of  laminar  combustion  based  on  con* 
sideration  of  reaction  rates  of  the  form  Ctxpi—TglT)  where  C»  1  is  a  nondimensional  preex> 
ponentiai  factor,  Tg  is  the  activation  temperature  and  T  is  the  temperature.  Typically  {Ta/T)»\  so 
that  the  exponential  factor  is  generally  small.  To  resolve  the  ambiguity  of  the  value  of  the  product  of 
the  two  factors,  one  large  and  one  small,  we  identify  a  temperature  Tg/lnC  so  that  the  reaction 
rate  becomes  exp((7'a/7t)(l~  in  this  form  we  see  the  product  and  thus  the  reaction  rate  is 

small  if  T <  Tk,  equal  to  unity  it  T  ^  Tk  and  large  if  7'>  Tk.  Thus  Tk  is  exposed  as  a  critical  tem¬ 
perature  while  the  reciprocal  of  the  temperature  ratio  Tg/Tk  is  an  appropriate  expansion  parameter. 

Since  its  development,  AEA  has  been  applied  to  a  variety  of  problems,  e.g.,  by  Law  [4],  Linan 
(5],  Khshnamurthy  et  al.  [6],  Peters  [7]  and  Kassoy  [8].  A  monograph  devoted  to  AEA  by  Buckmaster 
and  Ludford  (9]  is  near  publication.  In  a  recent  paper  Matalon  [10]  applies  AEA  to  the  description  of 
the  departure  from  chemical  equilibrium  in  the  gas  surrounding  a  carbon  particle  undergoing  oxidation. 
In  [2]  Libby  applies  this  method  to  the  burning  of  carbon  particles  in  order  to  establish  a  useful  ignition 
criterion  and  to  explain  qualiutively  the  abrupt  changes  in  particle  temperature  predicted  by  numerical 
solutions  to  the  complete  equations  of  thermal  response. 
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BACKGROUND  AND  ANALYSIS 

We  start  our  discussion  with  the  development  in  [2]  based  on  the  more  complete  analysis  in  [1];  these 
references  should  be  consulted  for  details.  It  is  sufficient  for  present  purposes  to  describe  the  basic 
situation.  We  consider  a  spherical  carbon  particle  of  initial  radius  r^o  suddenly  immersed  in  an  oxidiz¬ 
ing  ambient  with  temperature  and  with  a  mass  fraction  of  oxygen  denoted  Y]..  Nitrogen  treated  as 
an  inert  diluent  is  the  only  other  species  far  from  the  particle.  Only  the  direct,  heterogeneous  reaction 
Ch-  1/2  Oj— *CO  is  considered  so  that  the  gas  mixture  involves  oxygen,  carbon  monoxide  and  nitrogen. 
More  general  situations  including  those  involving  equilibrium  composition  in  the  gas  phase,  additional 
heterogeneous  reactions  and  additional  species  such  as  hydrogen  and  water  (cf.  [3])  can  be  treated  by 
the  methods  used  here,  although  increased  algebraic  complexity  is  involved.  The  temperature  T"*  is 
assumed  to  be  sufficiently  high  so  that  ignition  takes  place. 


The  Basic  Problem  and  its  As}'mptotic  Treatment 

The  physical  situation  described  here  leads  to  the  following  equations  which  are  written  in  a  form 
appropriate  for  asymptotic  analysis; 

fl'(T)- -|j- (1) 

(2) 

Af-eexp  j- i  l]  R  [(Yi.-h/i)  (3) 


A~  A(K)^  k[\-Cp  (l-e-^)-'j 
D  -  DiK)  -  (e-*  c,  (1  -  +  £] 


where 
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(  A,  M  J 

Equations  1-3  are  subject  to  the  initial  conditions  9(0)  -  (0)  - 1. 

The  principal  dependent  variables  in  these  equations  are:  9— the  particle  temperature 
divided  by  the  temperature  in  the  ambient;  K^Tphf, ,  the  instantaneous  radius  divided  by  its  initial 
value;  and  iO  0,  a  nondimensional  parameter  proportional  to  the  rate  of  mass  loss  of  the  particle. 
The  independent  variable  r  is  a  nondimensional  time.  The  other  quantities  in  Eqs.  1*3  are  thermo¬ 
chemical  parameters.  Ip,  A,,  A  and  the  various  /t's,  are  ccmstants  of  0(1)  related  to  the  thermo¬ 
dynamic  and  molecular  properties  of  the  solid  and  gas.  The  quantity  >  0  is  a  nondimensional 
parameter  related  to  the  radiative  properties  of  the  particle  and  the  ambient. 

Although  Eqs.  1-3  apply  to  a  wide  variety  of  conditions,  the  numerical  examples  in  the  ivesent 
study  involve  the  following  representative  values:  ^,<-0.186,  0.2,  A*4/3,  0.583,  £—2.18. 

More  significant  for  our  considerations  is  the  small  parameter  c  related  to  the  large  preexponential 
factor  k\  in  the  description  of  the  single  heterogeneous  reaction  accortUng  to 

ee'/*-£i  (4) 

where  £]  »  1  involves  the  product  of  a  constant  preexponential  factor  and  of  (prjx,/i7»)  Miere  p  is 
the  pressure  and  tj.  is  the  viscosity  of  the  gas  mixture  in  the  ambient.  Thus  increases  directly  with 
increases  in  pressure  and  particle  size  and  slowly  with  decreases  in  r«.  Note  that  conditions  of  practi¬ 
cal  interest  in  entrained  flows  meet  the  requirement  of  ^]  »  1. 

The  parameter  c  is  also  related  to  a  particular  value  of  the  nondimensional  temperature  denoted 
9t  which  is  the  critical  temperature  alluded  to  earlier  and  which  we  shall  find  is  the  ignition  tempera¬ 
ture;  the  relation  is: 


e/a  -  9*  -  0(1) 


(5) 


where  a=  TmUa,  the  ratio  of  the  ambient  and  activation  temperatures,  taken  to  be  fixed. 
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Several  comments  regarding  Eqs.  4  and  5  are  indicated.  For  a  specific  calculation  of  particle 
behavior  the  temperatures  T„  and  Tg  and  thus  the  parameter  a  are  known.  In  this  connection  we  note 
that  although  the  asymptotic  analysis  is  carried  out  for  »«>,  e— *0,  a  is  a  key  small  parameter 
determining  the  applicability  of  the  asymptotic  analysis  to  the  specific  calculation  in  question.  We  note 
that  an  alternative  asymptotic  analysis  could  involve  e,a— *0,  fixed  implying  that  both  k\  and  Tg 
increase  indefinitely.  Here  we  adopt  the  more  applied  point  of  view  and  imagine  that  we  are  interested 
in  the  effect  on  particle  behavior  of  changes  in  the  the  factor  In  a  specific  calculation  that 

factor  and  thus  k\  are  also  known  so  that  Eq.  4  determines  c  and  Eq.  5  determines  «4iich  may  thus 
from  our  adopted  point  of  view  be  considered  a  function  of  k\  and  a.  Increases  in  due  to  increases 
in  the  factor  (pr^/y),^)  lead  to  decreases  in  e  and  in  »* .  If  we  associate  with  an  ignition  tempera¬ 
ture,  we  conclude  that  the  temperature  at  which  vigorous  heterogeneous  reaction  is  initiated  decreases 
as  the  overall  reaction  rate  is  increased,  a  physically  appealing  conclusion.  However,  it  must  be  noted 
that  Eq.  4  calls  for  large  changes  in  to  achieve  relatively  small  changes  in  e.  To  illustrate  we  show 
in  Table  1  numerical  values  of  e,  ki  and  6*.  We  observe  that  a  halving  of  e  requires  an  increase  in  ki 
by  nearly  four  orders  of  magnitude  while  Eq.  5  implies  an  accompanying  reduction  of  by  only  one- 
half  'with  a  held  fixed. 

The  relations  involving  e,  k\  and  0^  influence  greatly  the  behavior  of  the  mass  loss  parameter  K 
via  Bj.  3.  A  consideration  of  the  range  of  possible  values  of  K  is  helpful  in  achieving  understanding  of 
the  asymptotic  analysis  and  its  applicability  to  a  specific  calculation.  To  that  end  we  rewrite  Eq.  3  as 


K 


c  exp 


R(Yx„  +  ii)  ie-'^  -e 


(6) 


where  the  special  value  of  K  denoted  Kj  is  defined  by  exp(-Arf)  —  +  and  is  the  diffusion- 

limited  rate  of  mass  loss,  the  rate  determined  by  the  ability  of  oxygen  to  diffuse  to  the  surface  of  the 
particle  independent  of  the  kinetics  of  the  heterogeneous  reaction.  For  the  values  of  the  thermochemi¬ 
cal  parameters  used  in  the  numerical  examples  0.140. 
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We  now  discuss  Eq.  6;  prior  to  ignition  when  0<9i,  (meaning  bounded  away  from  in  the  limit 
c— *0)  and  =  1,  Eq.  6  indicates  that  K  is  basically  an  exponentially  small  quantity  and  the  rate  of 
mass  loss  is  kinetically  controlled,  i.e.,  diffusion  of  oxygen  to  the  surface  of  the  particle  provides  no 
effective  inhibition  to  mass  loss.  Ignition  occurs  when  0;^— 9»0(c)  so  that  AT -0(c).  The 
identification  of  Bk  with  ignition  is  now  clear. 

Further  increases  in  B  above  Bk  lead  to  0(1)  values  of  K.  For  example,  wlien  iBk/B)  — 
1  — cln(l/c)+0(c)  and  /?— 0(1),  the  exponential  term  in  Eq.  6  is  large,  of  0(l/c)  and  Af-Od). 
However,  the  growth  of  K  is  limited  by  the  last  factor  on  the  right  side  of  Eq.  6,  i.e.,  by  AT— *A(j~,  so 
that  the  growth  in  the  magnitude  of  the  exponential  term  with  increasing  B  is  compensated  by  decreas¬ 
ing  values  of  the  last  factor.  When  K  —  Ka~,  Eq.  6  can  be  rewritten  in  the  following  illuminating 
form: 


A:/Arrfs=l-|eA/?exp||  1  - 


1-1 


(7) 


Equation  7  applies  provided  the  second  term  on  the  right  side  is  suitably  small  compared  to  unity.  This 
is  the  case  as  c— *0  provided  B  is  sufficiently  large  for  the  operative  order  of  magnitude  of  R.  For 
example,  if  /?  -0(1)  and 


^ _ 1 _  (0) 

Bk  1-6  In (/3(e)) 

where  /3»0(l/c),  e.g.,  /3— e“^  K^K^.  However,  as  complete  consumption  of  the  particle  is 
achieved  so  that  /?  -0(c),  then  Eq.  7  applies  and  the  rate  of  mass  loss  is  diffusion-limited  only  if 

_0_  _  1 

Bk  1  —  e  In  (yfc)) 


where  >>>0(1/6*),  e.g.,  y-c”*.  Compared  to  the  result  for  /? -0(1)  we  tee  that  the  temperature 
deviation  from  Bk  must  be  larger  when  /(  «  1  if  diffusion-limited  behavior  is  to  prevail.  It  shoul  d  be 
noted  that  when  0>0*  and  c-»0,  the  factor  cexp((l-(0*/0))/e)  is  expommtially  large  so  thjit  R 
must  be  exponentially  small  if  AT- 0(1)  <  A(,f.  In  contrast  we  observe  from  Eq.  7  that  for 
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/?“0(1)  K—^Ka  if  “ ein(y(«)/€/?)  where  y(«)  »  0(1)  and  R  have  an  algebraic  dependence 
on  €.  We  thus  conclude  from  Eq.  7  that  diffusion-limited  behavior  occurs  for  both  i?-  0(1)  and 
/?-0(l)  in  an  algebraic  sense  if  9-9^  -  0(<;in/ii(«))  provided  is  suitably  but  algebraically  large 
in  the  limit  e  — 0,  suitably  being  determined  by  the  requirement  that  the  second  term  in  Eq.  7  be  small 
compared  to  unity.  Thus  only  a  small  deviation  from  is  required  to  achieve  diffusion-limited  mass 
loss.  These  arguments  represent  a  generalization  of  the  9„-  sequence  introduced  by  Libby  (2]  and  used 
to  explain  qualitatively  the  particle  behavior  predicted  by  exact  numerical  solutions  to  the  describing 
equations. 

The  question  arise  as  to  the  applicability  of  these  results  related  to  c-*0  for  a  specific  calculation 
involving  a  finite  but  small  value  of  e.  From  a  consideration  of  Eqs.  3  and  7  we  cannot  be  assured  that 
for  a  particular  value  of  «  the  actual  histories  ^(t)  and  R{r)  respect  the  requirements  for  diffusion- 
limited  behavior  predicted  by  the  asymptotic  analysis  for  small  increments  in  9  above  9^ .  However,  we 
seek  within  the  asymptotic  analysis  itself  symptoms  of  failure,  i.e.,  failure  to  achieve  diffusion-limited 
behavior.  An  indication  of  such  failure  is  indicated  in  Table  2  where  we  show  the  dependence  of  the 
second  term  on  «  and  d  for  /?  - 1  with  the  thermochemical  values  used  in  the  numerical  examples. 

Note  that  for  our  specific  value  of  o,  namely  a -1/9,  the  limiting  value  of  e  for  ignition 
corresponds  to  ft*  =  1  and  is  thus  equal  to  0.111.  We  see  from  Table  2  that  the  range  of  values  of  e 
and  ft  yielding  diffusion-limited  behavior  in  a  specific  calculation  and  thus  consistency  with  the  asymp¬ 
totic  analysis  may  call  for  values  of  e  smaller  than  are  usually  found  to  be  acceptable  in  asymptotic  anal¬ 
yses. 

The  history  of  a  carbon  particle  suddenly  immersed  in  a  hot  oxidizing  ambient  as  described  earlier 
can  now  be  reexamined.  During  pre-ignition  ft  <  ft* ,  /?  -  0(1)  and  A"  <  exp  (-«~'(ft*/ft-  D),  i.e.,  K 
is  exponentially  small.  In  the  neighborhood  of  ft -ft*  abrupt  changes  take  place  in  ft  and  K,  changes 
associated  with  the  ignition  and  post-ignition  periods  which  are  found  to  have  durations  of  0(()  and 
0(«ln(l/«)).  Under  appropriate  circumstances  an  extensive  period  of  diffusion-limited  behavior  can 
ensue  with  ft  -  ft*  >  0(e  In  (l/e)).  Although  it  is  characteristic  of  all  of  these  periods  that  ft'  (r)  >  0, 
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eventually  the  right  side  of  Eq.  1  must  vanish.  Thereafter  a  complex  extinction  process  evolves  in 
which  9'  (t)  <  0  and  9  — •  1 . 

With  the  basic  situation  described  we  now  can  take  up  the  details  of  the  analysis.  In  [2]  the  pre- 
ignition  and  ignition  periods  are  studied  by  eliminating  r  between  Eqs.  1  and  2  so  that  the  differential 
equation  dR/dO  is  subjected  to  asymptotic  analysis.  In  this  approach  the  time  variation  is  found  in  the 
form  T  — t(9)  from  knowledge  of  R^Ris).  Here  we  find  analysis  of  the  particle  history  prior  to 
extinction  is  handled  more  conveniently  with  r  retained  as  the  independent  variable. 


Pre-lgnitlon  and  Ignition  Periods  . 

If  /?  —  0(1)  and  9(0)  <  9  <  9j^  so  that  AT  and  /?  - 1  are  exponentially  small,  *  1  and  Eq.  1  becomes 

9’(T)-3f,  (l-9)  +  0,(l-9*)  (9) 

which  describes  the  thermal  response  of  a  spatially  homogeneous,  inert  sphere  subject  to  conduaive 
and  radiative  heating.  A  supercritical  response,  i.e.,  one  involving  significant  chemical  reaction  and  par¬ 
ticle  consumption,  evolves  only  if  9'  (t)  >  0.  This  inequality  is  met  if  9*  <  1  so  that  we  esublish  one 
limitation  of  the  asymptotic  analysis.  For  a  given  value  of  Ki  and  thus  of  c  the  value  of  9/^  depends  on 
a;  we  thus  appreciate  its  important  role  in  determining  particle  behavior.  An  inverted  quadrature  of 
Eq.  9  gives  t*-t(9).  It  follows  that  the  duration  of  the  induction  period  given  by  t*  is  known  once  9* 
is  known  from  the  specified  values  of  e  and  a.  The  solution  of  Eq.  9  in  the  neighborhood  of  t  -t*  is 

9(t  <  t*)  *  9*  j  1  -1-  X^(t  -  t*)  -I-  •  •  •  I  (10) 

where  X*-3(f,  (1  -9*  )-f  (1  -9*  ))/9*  >  0  for  9*  <  1  has  a  known  value  assumed  to  be  0(1)  in  the 

limit  €—'0. 

Centered  about  t*  is  the  ignition  period  during  which  9  and  K  increase  and  R  decreases.  To 


describe  this  period  we  find  that  the  appropriate  scaling  transformations  are 
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e  =  (1  +  e  e,  +  •  •  •  ) 

R-l-e^Ri  (11) 

T  =•  T*  +  ef  . 

Then  Eo,  3  gives  Tioo  exp  (e, )  Tio„ «  1  to  a  first  approximation.  Wih  AiK)  and  DiK) 

appropriately  expanded  about  0  =  ^*  Eqs.  1  and  2  become 

0\  (t)  = 

^'2(t)  TiocC  *'  . 

The  following  solutions  which  match  properly  to  the  pre-ignition  solution  for  9  and  to  /?  *  1  are  readiiv 
found  to  be 


(t)  “  X*t 


i\2) 


ft  /-"i  y\oo  v2* 

R^ir)  « 

The  exponential  growth  of  K  limits  the  validity  of  the  expansion  in  Eq.  11.  In  particula’  whcr, 
€  exp  (X^t)  “=  0 (1),  then  A”*  0(1),  €^^2*  Ofe),  *=  €  In  (I/f)-!- 0(c)  and  cr  == 

€  In  (l/€)/X^+  0(€).  If  we  define  r,  *  X'^lnd/e),  then  it  follows  from  Eq.  12  that 

^(t  — •  T,)  —  1 1  —  6  In  (1/e)  +  X^  (t  “  T,)  -f  *1  (13) 

r,  ^  Tf,  -  T, 
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Post'IgnitloD  Period 

The  nonuniformities  in  the  solutions  for  the  ignition  period  imply  that  in  the  next  region 

jl  +  (elni)j  +  e«,  +  •  ■  •  ) 

^  1  -  e  ^,+  •  •  •  (14) 

Substitution  into  Eqs.  1-3  yields 


(15) 

It  should  be  noted  here  that  we  are  considering  a  particle  with  a  temperature  still  close  to  the  ignition 
temperature  and  with  a  radius  with  /i  ^  1.  However,  the  third  of  Eqs.  15  implies  that  significant 
mass  loss  is  encountered  because  0(1). 

As  long  as  the  term  within  brackets  on  the  right  side  of  the  first  of  Eqs.  15  is  positive,  the  tem¬ 
perature  variation  is  monotonically  increasing.  Accordingly,  it  worth  examining  the  variation  of  this 
factor  with  K.  In  Table  3  we  give  for  the  thermochemical  values  used  in  the  numerical  examples  the 
variation  of  the  quantities  A(/C)  and  D{K).  We  deduce  therefrom  that  the  factor  in  question  is  posi¬ 
tive  over  the  entire  range  of  K  provided  that  <  I  as  required  earlier. 

Equations  15  must  be  solved  numerically  subject  to  matching  to  Eq.  13  and  to  the  vanishing  of  k] 


as  r  %  -  eo;  thus  we  require 
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B]  (f  —  ~Oo)  ■=r  f 

A]  (f  —  -oo)  =r  0  . 

As  r— >oo,  0,  increases  in  an  unbounded  fashion  and  it  follows  that  as  a  consequence 

0',  (f  _  oo)  -  0.2  =  \^-A,9,  -I-  Ar  +  »,  (1  -  ) ) 

Aj  =  -A  (Ka)  >  0 
D^mD(Kd)  >  0 

so  that 

0(t  — •  oo)  -  0*  jl e  In  C-^) -I- e  (ft)2f -I- ATj) -I-  '■'j  • 

The  present  analysis  is  valid  provided  {u^-iOd)  >  0.  In  fact  for  specific  thermochemical  values  we  can 
determine  a  special  value  of  Bh  denoted  B^  which  makes  <o^  vanish.  Thus  we  have 

~AiiBi(  +  Dd  -f  (1  —  fl**)  “  0 

For  the  values  used  in  the  numerical  examples  — 1.234.  This  value  identifies  a  second  limit  on  the 
applicability  of  the  present  asymptotic  analysis  to  a  specific  calculation.  Physically,  Bl  is  the  tempera¬ 
ture  a  particle  with  /{  -  1  can  achieve  when  it  undergoes  diffusion-limited  chemical  reaction.  The 
analysis  is  self-consistent  provided  the  second  term  on  the  right  side  of  Eq.  7  with  B’^Bl  and  ^  1  is 

suitably  small  compared  to  unity.  This  requirement  defines  loosely  a  maximum  value  of  c;  for  exam¬ 
ple,  if  we  permit  the  second  term  to  be  0.1,  e  can  be  as  large  as  0.085.  If  we  require  that  term  to  be 
0.01,  c  can  be  as  large  as  0.07.  Accordingly,  when  we  compare  the  prediction  of  the  asymptotic  analysis 
with  exact  numerical  calculations,  we  take  c- 0.060  but  to  illustrate  the  nature  and  extent  of  the 
disagreement  which  results  when  the  asymptotic  analysis  is  overextended,  we  also  present  results  for 
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<>•0.090.  If  c  is  suitably  small,  we  conclude  that  once  a  supercritical  event  is  initiated,  particle  con¬ 
sumption  proceeds  with  vigor. 

The  interesting  nonuniformity  occurs  when  <r-  0(1).  Equations  14  then  imply  that  0(1) 
and  <01-  0(1).  Since  9\  becomes  unbounded,  it  follows  from  the  third  of  Eqs.  IS  that  K—'K^i  and 
that  the  diffusion  limit  is  achieved  while  R  —  0(1). 


Diffusion-Limited  Period 

The  nonuniformity  at  the  end  of  the  preceding  period  suggests  that  the  appropriate  variables  for  this 
succeeding  period  are 

9~9(s),  R-R(.s),  T-T*-l-^lnUj-t-^j,  “  1  “ 


€  fi  R 


(16) 

where  9>9ic.  It  follows  that  the  mass  loss  is  basically  constant  while  the  dependent  variables 

9(s)  and  R(s)  retain  their  primitive  form  implying  0(1)  changes  from  9^^  and  unity  respectively.  The 
equation  for  /{(s)  can  be  integrated  directly  to  yield 

)?-(1-25)''^  (17) 

where  we  impose  the  condition  /?  (s  —  0)  - 1.  Equation  1  becomes 

(->4rf9  +  Dd  +  «^,(l-2s)‘^^  (1-9^))  (18) 


Equation  18  must  be  solved  numerically  subject  to  the  initial  condition  9(5  -  0)  -  9*. 

We  obtain  immediately  from  Eq.  17  an  estimate  for  the  extinction  time  since  R-O  when  s  -  1/2; 
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T,  »  T*  +  (2  Ki)-'  (19) 

It  is  interesting  to  note  that  the  increment  added  to  the  induction  time  depends  only  on  the  concentra¬ 
tion  of  oxygen  in  the  ambient  and  on  the  molecular  weight  parameter  4-  This  implies  that  particle  life¬ 
times  in  pure  oxygen  should  be  one  quarter  those  in  air  provided  the  induction  period  is  relatively 
brief.  The  experimental  results  reproduced  in  [1]  for  the  lifetimes  of  a  wide  variety  of  carbon  contain¬ 
ing  particles  in  air  and  pure  oxygen  are  in  accord  with  this  estimate. 


Comparison  of  Exact  and  Asymptotic  Solutions 

The  combination  of  the  pre-ignition  period  plus  the  diffusion-liniited  period  obtained  from  the  solution 
of  Eq.  18  provides  the  main  features  of  particle  behavior  according  to  AEA.  At  this  juncture  it  is  thus 
appropriate  to  compare  exact  numerical  calculations  with  the  predictions  of  the  asymptotic  aiudysis.  As 
indicated  earlier  we  do  so  for  two  values  of  e:  0.060  and  0.090.  One  is  suitably  snudi  so  that  the  two 
limits  on  the  applicability  of  the  asymptotic  analysis  are  respected  and  the  second  too  large  for  the  con¬ 
dition  resulting  \n  K—'Kj  to  apply.  However,  both  values  of  c  with  our  fixed  value  of  a  yield  suitable 
values  of  namely  0.S40  and  0.810  respectively.  The  previously  given  values  of  the  thermochemical 
parameters  apply. 

Figures  1-3  give  the  variations  with  r  of  the  particle  temperature  parameter  9,  of  the  particle 
radius  in  terms  of  R  and  of  the  mass  loss  parameter  K  respectively.  For  the  larger  value  of  c,  i.e.,  for 
the  smaller  chemical  kinetic  rate,  ignition  is  delayed,  the  rate  of  mass  loss  is  smaller  than  the  diffusion- 
limited  rate  and  the  particle  lifetime  is  increased. 

Also  shown  on  these  figures  are  the  predictions  based  on  the  asymptotic  analysis  for  the  two  dom¬ 
inant  periods  of  particle  lifetime,  the  pre-ignition  and  diffusion-limited  periods.  First  we  call  attention 
to  the  response  of  the  inert  particle  as  given  by  the  solution  of  Eq.  9  and  the  two  values  of  the  ignition 
time  Ti,  associated  with  that  resptmse.  These  predictions  are  in  good  agreement  with  the  exact  solu¬ 
tions.  For  the  smaller  value  of  c  we  see  that  the  asymptotic  analysis  predicts  many  features  of  particle 
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behavior  with  an  accuracy  acceptable  for  most  purposes.  This  is  not  the  case  for  c<- 0.090.  Here  the 
difference  between  the  exact  and  diffusion-limited  rate  of  mass  loss  predicted  by  the  asymptotic  analysis 
leads  to  a  considerable  error  in  particle  lifetime.  It  is  worth  noting  that  for  e— 0.060  the  value  of  K 
given  by  the  exact  solution  is  slightly  less  than  reflecting  both  the  large  but  finite  magnitude  of  the 
exponential  factor  in  Eq.  3  and  the  small  but  nonzero  value  of  the  second  term  in  Eq.  7. 

It  will  be  found  convenient  in  the  analysis  of  the  extinctitm  period  to  obtain  solutions  in  the  form 
9-0(A)  and  it  is  therefore  of  interest  to  examine  the  details  of  the  transition  from  inert  particle 
behavior  (A  =0)  to  diffusion-limited  behavior  (A=s A/)  in  this  form.  We  do  so  for  both  the  exact 
and  asymptotic  solutions  in  Fig.  4  for  c— 0.060.  The  satisfactory  agreement  between  the  two  solutions 
is  noted. 

We  thus  conclude  from  these  results  that  AEA  provides  a  useful  approximate  description  of  parti¬ 
cle  behavior  provided  that  the  reaction  rate  paramei-er  Ai  is  sufficiently  large  for  the  expansion  parame¬ 
ter  e  to  be  appropriately  small  and  that  the  temperature  ratio  parameter  a  is  similarly  small. 


The  Extinction  Period 

Although  the  main  results  of  practical  utility  are  given  by  the  analysis  presented  to  this  juncture,  the 
elaborate  struaure  of  the  extinction  period  during  which  9—1  and  A— 0  is  interesting.  It  is  clear 
from  the  form  of  Eq.  18  that  a  singularity  develops  in  the  0-solution  «4»en  s  — ‘A”.  By  elementary 
means  we  find  that 


0 


^  0*  -  m 


(20) 


where 


m  =  3y/iil>r(9’*- ])/OAi- Ki)  . 


The  singularity  indicated  by  Eq.  20  is  weak  in  the  sense  that  0  approaches  a  finite  value  wliile  its 


Figure  4.  The  transition  from  inert  of  diffusion-limited  behavior  for  c- 0.060  (see  Fig.  1  for  legend). 


Burning  Gubon  Particles 


-20' 


derivative  is  infinite  and  defines  one  side  of  a  cusp-like  variation  of  0(r).  Values  of  practical  interest 
imply  that  m—  0(1)  >0  and  we  therefore  assume  this  inequality  prevails.  The  implication  is  that  as 
complete  consumption  of  the  particle  is  approached  its  temperature  increases. 

The  numerical  solution  corresponding  to  c- 0.060  and  shown  in  Fig.  1  indicates  the  nature  of  the 
0(1)  transition  made  by  9  from  0*  0*.  The  temperature  parameter  0*  — 1.7S  in  the  present  calcula¬ 

tions.  It  is  essential  that  0*  >  1  in  the  limit  e— ‘0  in  order  for  cusp-like  behavior  to  prevail.  We  note 
from  the  third  of  Eqs.  16  that  the  exponential  term  is  exponentially  small  as  c— *0  provided  0*  >  0|, .  It 
follows  that  R  must  likewise  be  small  before  the  mass  loss  rate  decreases  from  its  diffusion-limited 
value.  This  behavior  is  typical  of  systems  involving  high  activation  energies  where  the  reaction  once 
initiated  proceeds  vigorously  until  the  reactant  supply  (in  the  present  case  proportional  to  the  surface 
area  of  the  particle)  is  nearly  eliminated.  This  argument  is  consistent  with  our  earlier  discussion  of 
Eqs.  3  and  7. 

On  physical  grounds  we  expect  the  extinction  period  to  involve  a  reversal  of  the  sign  of  the  slope 
of  0(r)  from  the  positive  value  prevailing  at  the  end  of  the  diffusion-limited  period.  During  this  period 
0  declines  from  0*  to  I***  while  K  decreases  from  Kj  toward  zero.  We  find  that  this  period,  although 
brief,  involves  a  complex  structure  which  is  more  conveniently  analyzed  by  a  reformulation  of  the  basic 
equations  as  follows:  If  /I »  /?  (K,  0)  is  found  from  Eq.  3  and  used  to  convert  Eq.  2  into  the  equation 

we  obtain  the  single  equation  for  0-0(/r) 


dK 


o'  D 

^  ^  bK 
AT  +  JiM,. 


which  becomes  after  some  calculation 


iL 

dK 


3[.40-l-Z)-l-<6,^(l  -0^))  (l  +  ^(1  -exp(Ji:- JCrf))-' 
K^-  {lK9j€9^)(A9^  -0^)) 


(22) 
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If  an  appropriate  solution  to  Eq.  22  is  obtained,  Eq.  21  supplemented  by  R^R(K,9  (K)),  yields  Kir) 
and  hence  in  terms  of  AT  as  a  parameter  the  other  time  histories,  Rir)  and  9(t). 

We  show  in  Figs.  5  and  6  the  exact  numerical  solutions  in  terms  of  9(K)  and  R  (K)  for  € —0.060. 
As  time  increases  K  decreases  from  its  limiting  value  Kj- 0.140.  The  radius  parameter  R  decreases 
monotonically  toward  zero  while  the  temperature  parameter  9  exhibits  the  more  complex  behavior  sug¬ 
gested  earlier.  Our  task  is  to  provide  an  asymptotic  analysis  approximating  these  characteristics  as 
e-0. 

To  construct  a  description  of  the  first  portion  of  the  extinction  period,  denoted  Zone  1,  we  follow 
a  two-step  procedure;  first  Eqs.  18  and  20  suggest  the  forms 

9^9' (l-x(t)e+  ■■■  )  (23) 

K^  Ka(l-x(e)^+  ■  ■  ) 

It  follows  from  Eq.  3  that  if  x(«)®  «  0(«) 


R^xifVKiK  (24) 

The  importance  of  the  inequality  on  x  is  emphasized.  As  a  second  step  the  form  of  x(c)  is  found  by 
retaining  in  Eq.  1  the  second  terms  in  the  expansions  of  AiK)  and  DiK)  about  K  —  Kj  and  by  requir¬ 
ing  that  these  terms,  previously  neglected,  be  retained  as  c— »0.  As  a  result  we  find  a  distinguished 
limit  if 


X  -  (ACrfZ/ie)'*  exp 


(25) 


We  note  that  this  result  indicates  that  the  perturbations  in  Eq.  23  and  in  the  particle  size  given  by 
Eq.  24  are  exponentially  small  in  the  limit  c->0  provided  9'  is  suitably  large.  This  restriction  is 
respected  for  sufficiently  small  values  of  c,  e.g.,  for  c- 0.060.  The  different  behavior  near  extinction 


o 


(^ODd 


Figure  6.  The  radius  history  during  the  extinction  period  (see  Fig.  5  for  legend). 
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shown  in  Figs.  1*3  for  the  two  values  of  c  suggests  as  might  be  expected  that  the  asymptotic  analysis  is 
inapplicable  for  the  larger  value  near  extinction  as  well  as  earlier  in  the  particle  lifetime 

With  Eqs.  23  and  2S  substituted  into  £q.  22  we  obtain 


1. 

dK'^  K 


_3_ 

P* 


4>r  -  1) 

K}K^ 


(26) 


where  Ca  =  A  'iK  -  Kj)9'  +  D'iK  -  Kj)  -  - 1)/(1  -  and  a  s  3  are  positive 

constants.  To  match  the  solution  of  Eq.  26  to  the  behavior  of  Eq.  20  as  s  %  l/2~  requires 


5(Ar-0)  =  mKe'KiK) 


The  solution  of  Eq.  26  is 


“  6*  (fl  +  1)  $’{a-  1)  KiK 


(27) 


In  Eq.  27  a  complementary  solution  ^  is  suppressed  since  with  a  >  1,  the  required  behavior  as 
K—'O  prohibits  its  retention.  The  coefficient  of  K~^  equals  (m/9’K^)  so  that  matching  is  assured. 

The  0-solution  obtained  from  Eqs.  23  and  27  descnbes  the  approach  with  increases  in  time  toward 
0*  and  its  departure  therefrom  as  K  increases.  The  concomitant  decline  in  the  radius  is  obtained  from 
Eq.  24.  Consequently  the  solution  which  rounds  off  the  cusp-like  behavior  given  by  the  preceding 
period  takes  place  on  a  small  time  scale  defined  by  t  -  (t*  -l-  (l/2Kj)  +  X“^e  In  (1/e))  —  0  (x^(€)). 

We  show  on  Figs.  S  and  6  the  asymptotic  solutions  as  dashed  lines  within  the  range  of  validity  of 
K  and  extend  them  beyond  that  range  as  dot-dashed  lines.  When  this  procedure  is  followed  for  the 
solutions  given  by  Eqs.  27  and  24,  we  find  that  within  plotting  accuracy  the  asymptotic  and  numerical 
solutions  carmot  be  distinguished  within  the  range  of  validity  of  the  former.  As  a  consequence  only  the 
extensions  of  these  solutions  for  Zone  1  appear  on  Figs.  5  and  6. 
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To  find  the  solution  in  the  layer  denoted  Zone  2  let 


(28) 


K  -^Ki  (I  -  e  it  +  •  •  • ) 


As  a  consequence  of  these  forms 


R  =•  x^exp  (fl*  f/fl' )/«  Kj  k 


which  indicates  a  further  dramatic  reduction  in  particle  radius.  Then  Eq.  22  becomes 


dk 


_ ag  -  (3  Cg/B')  k _ 

it(l  -  (0J9')  {ag~  OQl9')k\ 


(29) 


which  is  to  be  solved  subject  to  the  matching  condition  g(k—^0)-*0. 

For  increasing  values  of  k  corresponding  to  increasing  time  the  asymptotic  form  of  Eq.  29  is 


dg  ^  9' 
dk  9^k 


which  has  the  solution 


^ (k  —  oo )  =r  In  k  +  .  (30) 

It  follows  from  Eqs.  28  and  30  that  the  solution  in  this  zone  is  singular  when  c/c^O  (1)  which 
corresponds  to  eg  -  0  (c  In  (1/c))  and  R  -0  (xVc)-  Significantly,  the  nonunifoimity  with  respect  to 
K  suggests  that  the  mass  loss  declines  significantly  from  Ka  in  the  next  zone. 

Numerical  solutions  to  Eq.  29  and  to  the  related  solution  for  R^‘R(K)  are  added  to  Figs.  S  and  6 
to  permit  comparison  with  the  exact  n  jmerical  results.  Again  the  asymptotic  solutions  are  continued  by 


Burning  Carbon  Particles 


-26- 


dot-dash  lines  beyond  an  appropriate  value  of  K  to  indicate  their  range  of  validity.  Satisfactory  agree¬ 
ment  between  the  two  solutions  is  achieved. 

We  find  that  in  the  next  zone.  Zone  3,  in  which  K  declines  rapidly  we  must  take 


1 -f-e  In  (-)-«/((/:) -I- 
9k  « 


(31) 


The  equation  for  h{K)  is 


dh  _ 

9i_ 

1 

dK 

9k 

l-exp(A- Arf)  ^ 

which  is  to  be  solved  subject  to  the  matching  condition 


(32) 


The  Ovi)  constant  of  integration  gg  is  obtained  from  the  #r-solution  (Eq.  29).  The  solution  of  Eq.  32  is 


A  -  (in  U  -  In  exp  ((/frf  -  .«:)  -  ll  1 

6»*  l  l  t  I  J 


(33) 


When  comparison  between  the  asymptotic  and  numerical  solutions  is  considered,  we  encounter  a 
practical  limitation  of  the  former  solutions.  For  the  values  of  the  parameters  used  for  the  numerical 
solutions  we  find  that  the  second  term  on  the  right  side  of  Eq.  31  is  0.54,  indicating  that  the  value  of 
c- 0.060  overstrains  the  asymptotic  analysis.  Accordingly,  we  cease  making  comparison  of  the  two 
solutions  but  for  completeness  carry  through  the  asymptotic  analysis  with  the  expectation  that  agree¬ 
ment  with  numerical  solutions  would  be  achieved  for  suitably  small  values  of  c.  For  example,  with  all 
other  parameters  fixed  at  their  present  values  a  value  of  c— 0.022  would  make  the  second  term  in 
Eq.  31  equal  to  0.1,  an  appropriately  small  value. 

When  Eq.  33  is  substituted  into  Eq.  32,  we  find  that  a  nonuniformity  occurs  when 
(9’l9k)t  In  (1/A)*-  0(1)  implying  that  the  mass  loss  rate  becomes  vanishingly  small  at  the  end  of  this 


Burning  Carbon  Particles 


-27- 


zone.  As  a  consequence  the  temperature  decreases  further  from  9*.  Physically,  the  chemical  reaction 
weakens,  permitting  convective  cooling  to  reduce  the  temperature  of  the  dying  particle  so  that  9-*l. 
Also  in  this  zone  Eq.  3  indicates  that  R  remains  nearly  constant  at  a  small  magnitude  of  0(xV<). 

The  mqjor  decline  in  9  toward  unity  occurs  in  Zone  4  in  which  it  is  found  convenient  to  scale  K 
by  the  transformation 


K  - 


9kP 

9't 


(34) 


which  defines  p  as  a  new  0(1)  independent  variable.  This  exponential  transformation  is  suggested  by 
the  nonuniform  behavior  of  the  solution  in  the  preceding  zone.  Then  to  all  algebraic  orders  in  e  the  9- 
equation  is  simply 


^  _  ii 

dp"  9’ 


(35) 


9(p  — '  0)  -  9* 


where  the  latter  is  the  appropriate  matching  condition.  The  solution  of  Eq.  35  is 


9- 


1  +p 


0  <  p  <  (9*  -  1) 


where  the  limit  on  p  is  required  to  prevent  the  temperature  from  dropping  below  unity.  During  this 
cooling  phase  the  exponentially  small  radius  is  basically  constant  and  given  by 


R^R.~ 


1 


exp 


-F 


(36) 


These  results  can  be  used  in  an  appropriate  form  of  Eq.  1  to  show  the  time-dependence  of  the  tempera¬ 
ture  decline;  we  find  that 


9  -  1  +  (9*  -  1)  exp 


(t-t*) 


(37) 


-28- 


Buming  Carbon  Particles 

where  T- .  r .  +  a/2X.)  +  x'^c  In  d/*)-  'H'e  small  variation  from  the  extinction  time  should^  not^, 

D  -^o  an<t  fl— I  This  xone  provides  little  additional 
Further  decay  of  a  minute  nature  occurs  as  R 

irfom^ion  on  botavior  «.<!  is  no.  consider..  f»r«»r. 


CONCLUDING  REMARKS 


The  method  of  activation  energy  asymptotics  is  applied  to  the  burning  of  carbon  particles  in  a  hot  oxid- 
ing  ambient.  To  simplify  the  analysis  gas  phase  reactions  are  assumed  negligible  and  only  the  direct 
heterogeneous  reaction  between  carbon  and  oxygen  is  taken  into  account.  It  is  shown  that  a  complete 
description  of  particle  history  is  largely  given  by  an  induction  period  during  which  the  particle  is  heated 
to  an  ignition  temperature  and  by  a  subsequent  period  during  which  the  rate  of  mass  loss  is  diffusion- 
limited.  Such  a  picture  is  not  dissimilar  to  earlier  approximate  calculations  of  particle  behavior  based  on 
ad  hoc  assumptions.  However,  the  present  analysis  provides  clear  indications  of  the  conditions  under 
which  the  analysis  applies.  In  fact  a  significant  portion  of  our  discussion  relates  to  the  applicability  to  a 
specific  calculation  involving  c  small  compared  to  unity  but  nonzero  of  the  limiting  behavior  arising 
when  the  small  parameter  c-*0.  The  asymptotic  analysis  provides  estimates  for  quantities  of  applied 
interest;  conditions  under  which  ignition  occurs,  conditions  for  diffusion-limited  mass  loss  and  the 
times  for  ignition  and  complete  consumption. 

The  results  of  two  representative  numerical  examples,  one  respecting  the  limits  of  applicability  of 
the  asymptotic  analysis  and  a  second  overextending  that  analysis,  are  compared  with  the  predictions  of 
AEA.  Good  agreement  is  found  for  the  former  while  as  expected  significant  disagreement  is  found  for 
the  latter. 

During  a  brief  extinction  period  all  of  the  quantities  describing  particle  behavior  undergo 
significant  variations.  As  a  consequence  the  asymptotic  solutions  for  this  period  are  complex  and 
impose  significant  limitations  on  the  acceptable  size  of  the  expansion  parameter.  In  fact  only  two  of  the 
four  zones  found  to  be  required  are  compared  to  an  exact  numerical  solutions. 

Extension  of  the  present  analysis  to  include  gas  phase  reactions  and  more  complex  composition  in 
the  surrounding  ambient  can  be  carried  out  at  the  expense  of  increased  algebraic  complexity.  However, 
the  main  features  of  the  analysis  of  applied  interest,  namely  the  description  of  the  induction  and  a 
diffusion-limited  periods,  may  involve  straight-forward  considerations. 
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nomenclature 

a,Cj,m  parameters  arising  in  the  extinction  period 
A  (K),  D(K)  functions  of  mass  loss  parameter 

nondimensional  parameter  related  to  the  coefficient  of  specific  heat  at  constant  pressure 

E  thermochemical  parameter 

g,  h  perturbation  functions  arising  in  extinction  period 

K  mass  loss  parameter 

Kj  diffusion-limited  value  of  K 

k\  nondimensional  reaction  rate  parameter 

k  perturbed  mass-loss  parameter 

r,  instantaneous  particle  radius 

/■po  initial  particle  radius 

R  nondimensional  radius  parameter 

5  transformed  time  variable 

T  temperature 

Ta  activation  temperature 

T„  temperature  in  surrounding  ambient 

Y\oo  oxygen  concentration  in  surrounding  ambient 

Greek  Symbols 

a  temperature  ratio,  Too/ Ti  ' 

A,  A,  thermochemical  parameters 
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9  nondimensional  temperature  parameter 

nondimensional  temperature  parameters  arising  in  AEA 
e  expansion  parameter 

/i,  fi.,  m'  molecular  weight  parameters 

radiative  parameter 
parameters  in  AEA 
r  nondimensional  time 

f  transformed  time  variable 
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